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Recently, Girolami and Adesso have demonstrated that the calculation of quantum discord for two-qubit case 
can be viewed as to solve a pair of transcendental equation (Phys. Rev. A, 83, 052108(2011)). In present work, 
we introduce the generalized Choi-Jamiolkowski isomorphism and apply it as a convenient tool for constructing 
transcendental equations. For the general two-qubit case, we show that the transcendental equations always 
have a finite set of universal solutions, this result can be viewed as a generalization of the one get by Ali, Rau, 
and Alber (Phys. Rev. A, 81, 042105 (2010)). For a subclass of X state, we find the analytical solutions by 
solving the transcendental equations. 
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I. INTRODUCTION 

How to quantify and characterize the nature of correlations 
in a quantum state, besides the fundamental scientific interest, 
has a crucial applicative importance in the field of quantum 
information processing [1]. For a bipartite quantum state, it 
is known that both the classical and quantum correlations are 
contained in it. Beyond the entanglement, quantum discord 
was introduced as a more general measure of quantum cor¬ 
relation [2, 3], and was regarded as a resource for quantum 
computation [4], quantum state merging [5, 6]. Quantum dis¬ 
cord has attracted much attention recently [4—15], and has also 
been generalized to continuous-variable systems for Gaussian 
states [16, 17] and non-Gaussian states [18]. 

Quantum discord is very hard to calculate even for two- 
qubit states because of the minimization over all possible mea¬ 
surements. For an important class of two-qubit states, the 
so-called X states, Ali, Rau, and Alber (ARA) proposed an 
algorithm to calculate the quantum discord with minimization 
taken over only a few simple cases [9] . However, a counterex¬ 
ample for the ARA algorithm was given by Lu et. al. [10], 
where the authors proved that, for the entire class os X states, 
the optimization procedure involved in the classical correla¬ 
tion should be state dependent. For the real X states, Chen 
et. al. have identified a class of states, where quantum dis¬ 
cord can be evaluated analytically without any minimization, 
and hence the ARA algorithm is valid. Meanwhile, they also 
identified a family of states for which the ARA algorithm fails 
[19]. 

The ARA algorithm involved a minimization procedure 
with four constrained parameters. However, Girolami and 
Adesso have shown that two free parameters, the polar and az¬ 
imuthal angles usually used to describe an arbitrary unit Bloch 
vector, are already sufficient. With the two angles, one may 
obtain two partial derivatives for the conditional entropy, and 
by setting the two partial derivatives to be zero, the minimiza¬ 
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tion procedure can be simplified as to find the solutions of 
a pair of transcendental equations [20]. Usually, one should 
firstly give all the possible solutions, which are series of val¬ 
ues of the two angles, and then select out the optimal setting 
where the conditional entropy takes the minimal value. 

Although the transcendental equations is direct and reliable, 
it has been argued that, for general case, one cannot solve 
the problem analytically since these equations involves log¬ 
arithms of nonlinear quantities [20]. In preset work, we shall 
give some further discussion about this problem. First, we 
introduce the generalized Choi-Jamiolkowski isomorphism as 
a convenient tool to construct the transcendental equations. 
Then, for the general two-qubit case, we demonstrate that 
the transcendental equations have a set of universal solutions 
which have been discovered by the ARA algorithm. Finally, 
for a subclass of the X state, we give the analytical solution 
by solving the transcendental equations. 

The content of present work is organized as follows. In 
Sec. IT, we give a brief review of the quantum discord. In 
Sec. III., we introduce the general Choi-Jamiolkowski isomor¬ 
phism. In Sec. IV, a detail introduction of the Bloch vector 
transformation is discussed. In Sec. V, we give a classifica¬ 
tion of the solutions for the partial equation of the classical 
mutual information. In Sec. VI, several examples are given 
there. Finally, we end our work with a short conclusion. 

II. THE QUANTUM DISCORD 

The correlations for a bipartite state can be quantified by 
the quantum mutual information. For a given density matrix 
pab ^ bipartite system C) the quantum mutual infor¬ 
mation is defined as 

I = 5(p“) + (1) 

where S{p) = —Tr(plog 2 p) is the von Neumann entropy, 
and p“ (p^) denotes the reduced density matrix of subsystem 
Ra xhe quantum mutual information can be expressed 

as the sum of two part. 


( 2 ) 
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with C(p°‘^) the classic correlation and Q{p°'^) the quantum 
discord [2, 3]. To quantify the quantum discord, Olliver and 
Zurek [2] has suggested the use of von Neumann type mea¬ 
surements: with Ili the one-dimensional projective 

operators. After the measurement on subsystem H^, a density 
operator pj associated with the outcome j is 

Pj = —(Id 0 n,)(p“'')(iD 0 n,), (3) 

Pj 

with Pj the probability for the j-th outcome. Use 5(^111^) = 
'YhPj^iPj) to denote the quantum conditional entropy, and 
the corresponding quantum mutual information reads 

i'(p“'|n,) = 5(p“)-5(p|n,). (4) 

Then, the classical correlation is 

C(p“'’):=sup{n,}l'(p“'|n,), (5) 

and the quantum discord is defined as 

Q(p“'') (6) 

III. THE SYSTEM-ANCILLA-ENVIRONMENT PICTURE 


well the entanglement between the two system is preserved 
by the quantum process e [22]. In the following, we shall 
show that this process is reversible: If the reduced matrix of 
is full-rank, it can always be described by the correspond¬ 
ing e and |<I>). To prove this, we should first notice that a 
bounded operator in Hd is always related to a vector in a ex¬ 
tended Hilbert space H®^. Denote A to be a bounded operator 
on the £)-dimensional Hilbert space H^, with = (i|A|j) 
the matrix elements, and an isomorphism between A and a 
-dimensional vector | A)) can be 

D 

|A)) = VdA(^1d\S+) = Y, (8) 

where |S'+) is the maximally entangled state in H®^, |S'-|-) = 

SfcLi \kk)/\/D with \ij) = |i) 0 \j). This isomorphism 
offers a one-to-one map between an operator and its vector 
form. Suppose that A , B, and p are three arbitrary bounded 
operators on H^, and then 

Tr(A+B) = {{A\B)),\ApB)) = A 0 B^\p)), (9) 

with the transpose of B. 

Now, consider a x density matrix. 


The Choi-Jamiolkowski isomorphism is a useful connec¬ 
tion between quantum channel and a bipartite state [21], say, 
= e 0 l£)(|S'+)(S'+|), with |S'+) the maximally entan¬ 
gled state of the bipartite system. Our work is motivated by 
such a simple idea: We first express a density operator 
with a quantum channel e, and then the analytic expression of 
the quantum discord for the D ® D system can be simplified 
since only a Z9-dimensional quantum process e is involved. It 
should be noticed that the isomorphism above is only avail¬ 
able for the cases when the reduced density matrix is a 
completely mixture, = 1d/D. With careful analysis, we 
find that the isomorphism above can take a general form as 
the maximally entangled state is substituted by a general en¬ 
tangled state, and then, the density matrix with a full-rank 
reduced matrix p^ can always be expressed with a quantum 
channel and an entangled state. Meanwhile, the Bloch vector 
transformation can be applied to describe the quantum opera¬ 
tion for the qubit case. Therefore, the derivation of the quan¬ 
tum discord is closely related to the property of the quantum 
channel. 

To study a quantum channel £ of a Z9-dimensional system 
iJ“, it is convenient to introduce an ancilla system with 
an equal dimension. Prepare a pure entangled state |<I>) as 
the initial state of the bipartite system 0 H^, and since 
the system is subjected to a interaction described by the 
trace-preserving quantum operation e with the environment, 
the final state is 


p“«' = e0l^(|$)($|). (7) 

From it, we can obtain a lot of information about the quantum 
channel. For example, the Schmacher’s channel fidelity is de¬ 
fined as F = (4>|p“^|<I>), which provides a measure of how 


^ ( 10 ) 

m—1 

where !'!'„) are the normalized eigenvectors of the bipartite 

density operator (T'm|4'n) = Smn, and Am are the cor- 

_^ £^2 

responding eigenvalues with Am = 1- Due to the iso¬ 
morphism I'km) = |rm)). the density matrix can be also 
expressed as 

p“' = ^Am|rm))((rm|, 

m 

and the transpose of the reduced density matrix p^ can be de¬ 
rived as 


(p'')T = ^AmrJ„rm. (II) 

m 


A simple proof of Eq. (11) is as following: With the equations 
in Eq. (9) and Eq. (10), one can obtain ((Fml^*) = (*|rLl^)’ 
and (fcjjFm)) = {k\Tm\j)- Plugging these results into the 
definition of the partial trace operation, we have 


mpY\j) = (jI(p^)i*) 

D 

k—1 m 

m k 


m 


□ 
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For the two-qubit case, D = 2, and one can further assume 
det(p^)"^) ^ 0, and define 

y^=t/F!C/t, (12) 

where 



with t/ a 2 X 2 unitary transformation. Furthermore, one can 
introduce a set of Kraus operators {Em}m=i for rh® quantum 
process e as follows; 


with pj the probabilities 

Pi = -(1 -I- cos0cos7),P2 ~ 2 ^^ ~ COS0COS7), (19) 
and \(f>j){j = 1, 2) are a pair of pure states defined as 

\(t>i) = ;^(cos|cos^|0)-t-sin|sin^e"*‘^|l)), 

\h) = ^(-cos^sin^lO)-f sin^cos|e"*'^|l)). 

2 2 2 2 

Finally, one can obtain 


Em — Ej^Em — I2, 

m 

and now, the density matrix p°‘^ can be rewritten as 

= Y.Em® i 2 m^)){m^\Ei 012 

m 

= Y,Em®u*mm\Ei®{u*)K 


(13) 


2 

^a'b' ^ '^p.p. 0 = e{\cj)j){(l)j\). ( 20 ) 

f=i 

Meanwhile, it is easy to check that Pj\4’j){4>j\ = ^rid 
therefore 

2 

=e{p^) = p^. ( 21 ) 

i=i 


It is obvious that a new basis for can be defined as 
|0) = (C/*)t|0),|l) = (C/*)t|l), 


and the relation in Eq. (7) with the entangled state 

1$) = \n)) = cos^lOO) -f sin||ll). 

can be obtained. In this picture, the two reduced density ma¬ 
trices are 


P 


b 


COS^ ^ 0 

0 sin^ ^ 


= e(P^)- 


(14) 


It should be mentioned that has the same determinant 

as p^, and furthermore, our method above can be easily gen¬ 
eralized for the cases with the arbitrary dimension D. 

In the following, we shall focus on the situation where a von 
Neumann measurement is performed on subsystem Two 
free parameters, 9 and </), can be used for the measurement 
operators Hi = \'ij}i){il}i\{i = 1, 2), where 


9 6 

IV'i) = cos-|0)-f sin-exp(z(/))|l), (15) 

6 6 

IV' 2 ) = -sin-|0)-f cos-exp(i(^)|l). (16) 


From Eq. (20), we see that the classic information I' is a 
function of the free parameters 0 and (f>, 

2 

I'{9,^)=S{p^)-'£p,S{p,), (22) 

2 = 1 

and the quantum discord can be accessed if the maximum 
value of X'{9, (j)) has been decided 

Q = I-Maxe.0l'(0,(/)). (23) 

IV. THE BLOCH VECTOR TRANSFORMATION 


In order to obtain the analytic expression of the quantum 
discord for a two-qubit state, we shall at first give a general ex¬ 
pression of the conditional entropy Ylj=i Pji^iPj)- Th® Bloch 
representation is very useful for the single-qubit state, and the 
state p can be written as p = | (I 2 + r ■ (?) with r is a three 
component real vector and a = {(Jx,cry,az)- Meanwhile, it 
turns out that an arbitrary trace-preserving quantum operation 
is equivalent to a map such that 

r' —>• r = pr-\- c, (24) 


After the measurement, the final state p“ ^ can be formally 
expressed as p“ = e 0 l 2 (p), where 

2 

P = ^i2 0nj|o))((o|i2 0n,. (17) 

2=1 


By some algebra, we find that p is a mixture of product state 


2 

2=1 


with p a 3 X 3 real matrix, c a constant vector, and e(p) = 
i (I 2 + r' • (f ). This is an affine map, mapping the Bloch sphere 
into itself [1], and can be explicitly expressed as 

( 9xx Vxy 9xz \ f ^ ^ \ \ 

Pyx Vyy Vyz 1 | I’y 1 + | Cy 1 ) (25) 

Vzx tlzy r]zz J \rz J V Cz / 



with the coefficients defined as 


Vij = ^Tr[crje(crO],Cfc = ]^Er[ake{l2)\. 


(18) 


(26) 
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Here, we have used 


= ^(l 2 +S-ct), \(j)2){(t)2\ = ]^{l2+t- a), (27) 
and the two unit vectors s and t 


Sz 

ty 


sin 7 sin 0 cos (j) sin 7 sin 6 sin (j) 

1 + cos 7 cos 9 ’ ^ 1 + cos 7 cos 9 ’ 

cos 7 + cos 9 ^ — sin 7 sin 9 cos (p 

1 + COS 7 COS 0 ’ ^ 1 —cos 7 cos 0 

— sin 7 sin 9 sin cp ^ cos 7 — cos 0 
1 —cos 7 cos 0 ’ ^ 1 —cos 7 cos 0 


(28) 


With the following two vectors, 

s' = rjs + c,P = r]t +c (29) 


one may have 


Pi = ^{h + s' ■a),p2 = ^{h+P ■ ct), 


(30) 


For simplicity, s'{9,p) and t'{9,(p) are used to denote the 
purity of the density matrix pi and p 2 respectively, and 

+ + t'{9,p) = n 

It is easy to note that there exists a symmetry between these 
two functions: Under the transformation 


0 —>■ TT — 0, ^ —>■(/) + TT (31) 

these two functions are interchanged 

s'{9,(P)^t'{9,(P). (32) 

This result comes from the fact that r (tt — 0, ^ + tt) = s(0, (p), 
which can be seen from Eq. (28). 


V. CLASSIFICATION OF THE SOLUTIONS 


With the Bloch vector introduced in above section, one may 
get a general expression of the classic information, 

I'{9,P) = (33) 

with H 2 {p) the binary entropy defined as 7 ^ 2 (p) = 

-plog 2 p-(l-p)log 2 (l-p).FromEq. (19), ^ 
Therefore, we can obtain 


dl' 

d(p 

dl' 


d(p 


{Pl log 2 


1 + s' dt' 

— ) + ^fclofe 


1 + t' ^ 
1-f'^ 






■)) 


(34) 


1 + s' dt' 

7) + ^(P2log2 


1 - s' 


1 + f' ^ 
1-t'' 


In the following, we shall show that the following two types 
of solutions are universal: 

{A)The symmetric solution: In this case, one of the solu¬ 
tions happens with the setting 



with p is constrained by 


s'{9,p) = s'{9,tt + P), 
ds'{TT/2,p) _ dt'{TT/2,p) 

dp dp 


(36) 

(37) 

(38) 


Following the discussions about the symmetry between 
s'( 0 , p) and t'( 0 , p), there should be 


s'{ttI 2,P) = t'{Trl2,p), 

ds'{9,p)^ _ 5f'(0,^)| 

--le= 77/2 - -le= 7 r/ 2 - 

Jointing the above results with 

pi(0 = tt/2) = i,p 2(0 = 7 r/ 2 ) = 


(39) 

(40) 


(41) 


one can conclude that the setting in Eq. (36) is one of the 
solutions. 

(B)r/!e asymmetric solution: Another solution of the partial 
equation exists with the setting 


9 = 9,P = P, (42) 

with 0 and p the solution of the equations below. 


ds'{9,p) ^ _ dt'(9, p ) , _ 

50 50 

4>) I _ dt'{9,p) 

dp dp 

--^1 --0 


(43) 

(44) 

(45) 


Except the special case where pi = P 2 = 1/2, 0 = 0 is the 
only possible solution for the equations above since dpi/d9 cx 
sin0, {i = 1, 2). By jointing it with the symmetric solution, 
0 = 7 r/ 2 , the main result in [9], which states that the polar 
angle 0 may take the value 0 or 7 r /2 for the X states, is also 
suitable for the general two-qubit case. 

{C)The state-dependent solution: For some given states, the 
two transcendental equations in Eq. (35) may have other type 
solutions beside the universal one given above. We shall give 
an example in the next section. 


VI. EXAMPLES 

(A) The X state. This type of density matrix has been 
widely discussed in previous works [7, 9, 10], 


As a necessary condition, the maximum value may happen 
with 


dl'/dp = 0 ,51750 = 0. (35) 


^ah 


/ Pit 0 0 pi4 \ 

0 P22 P23 0 I 

0 P32 P33 0 I 

V P41 0 0 P44 / 


(46) 
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By some simple algebra, we may see that the map in 
Eq. (25) now take the form 


rn 1 

1 ^XX 

Ixy 

0 


Vyx 

^yy 

0 

\r'z 

\ 0 

0 

qzz 




the equation in (55) has no other solutions beside 0 = 0 and 
s' = t'. From Eq. (56), there should be 


dH{x) 


1 


dx + kx'^ 

with the expanding formula 


[h 


1 


1 — X 


(58) 


1 — 


Here, we focus on the case when cos 7 = 0, which means 
the pure state |$) in Eq. (7) is the maximally entangled state. 
Under this condition, the probability for each final state takes 
the same value, pi = P 2 = 1 /2. The transcendental equations 
are reduced as 


(1 + 

1 — 


= 2x(l + 

n—1 


2'/7' -\- 1 


= l + (l + fc)^x2^ 

n—1 


0 

0 


d4> 


(log2 


1 + I 


1 + f', 

l-f' 


ds' 


1 + s' dt' 


l + t' 


l-t 


l)- 


(48) 

(49) 


With the vector transformation, we shall get 

s'(6»,(/)) = {sm^6»/((/)) + [c^+P2 ^cos 6']^}^, (50) 

= {sin^^'/(</') + [c^-Pzzcose]^}^, (51) 

f{4’) = iVxx cos 4> + -q^y sin (j)f + {qyx cos + qyy sin (j)f' 


and the condition in Eq. (57), we see that dH (x) jdx is non¬ 
zero in the parameter range 0 < x < 1. Therefore, H{s') = 
Hit') can only happen with s' = t' . If Cz ^ 0, s' = t' has 
the unique setting 0 = 7r/2. If = 0, s' = f' can hold for 
an arbitrary 0, while from Eqs. (50-51), we find the optimal 
setting is either 0 = Oor0 = 7 r/2 . Based on these analysis 
above, we conclude that the universal solutions are sufficient 
for the cases above. 

Among all the W-type states, the Bell diagonal state is one 
of the most interesting cases, and in the parameterized state 
model here, it corresponds to the situation 


Note that s' and t' depend on (j) in the same way. There¬ 
fore, the optimal setting for ip should be decided by equation, 
df{(j))/d(j) = 0, can be easily solved. By introducing a set of 


parameters. 




ril = 

max/((/)), a = 
<P 


(52) 

c = 

2 2 7 

Vzz = 

C 

(53) 

b"^ — ca 


we can express s' and t' with a simple form, 

s' = {a-|-26cos0-|-ccos^ ,t' = {a—2bcos0+ccos^ 0}^. 


From it, we get the derivatives. 


ds^ /- 

l-^l = sm6»A/62 - caG(s'), 

11^1 = sm6»v/62 - caGit'), 

G(x) = (54) 

X 

Note that two derivatives can not be positive at the same time, 
we can rewrite Eq. (49) as 


0 

Hix) 


sm0\/b'^ — ( 


-{His')-Hit')), 


21 n2 

v/l + kx'^ , 1 + a; „ , 

-In--,0 < X < 1 

X 1 — X 


(55) 

(56) 


qxx 0 ^ \ ^ \ 

0 qyy 0 r, (59) 

0 0 qzz J \rz J 


Now, the parameter k takes the value k = —1/ql < —1. The 
symmetric solution should be s'(7r/2, = f'(7r/2, cp) = q± 

while the asymmetric solution has a compact form s'i0,(p) = 
t'i0, cp) = \qzz\- Finally, which kind of solution, the symmet¬ 
ric one or the asymmetric one, should be viewed as the classic 
correlation C in Eq. (5), is decided by the actual values of q^x, 
qyy,and qzz- Formally, it an be expressed as 

C = l_772(i^|2Pi)^ (60) 

with qopt — I j 1 ^ 2.2 |}- 

(B) In Ref. [10], a simple density matrix is given as 

/ 0.0783 0 0 0 \ 

afc_( 0 0.1250 0.1000 0 I 

^ ~ 0 0.1000 0.1250 0 ■ 

Vo 0 0 0.6170 / 


With numerical calculation, we find the transcendental equa¬ 
tions in Eq. (35) have three solutions, 0 = 0, 0 = -k 12 and 
0 K, O.ISStt. Among all these possible settings, 0 k, 0.1557r 
is the optimal one. With this simple example, we show that 
the universal solutions are not always the optimal one. 


VII. CONCLUSIONS 


Here, we shall show that: In the parameter range 

2 

k > —- or k < —1, (57) 

o 


Our present work has followed the original definition of the 
quantum discord in Ref. [2], where the von Neumann projec¬ 
tive measurement is performed. This measurement can also 
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be generalized to the more general positive operator-valued 
measurement (POVM) [3]. Furthermore, the concept of the 
quantum discord itself has been developed in recent years. For 
examples, the relative entropy quantum discord [12], the geo¬ 
metric quantum discord [13, 14] and their relations to the orig¬ 
inal definition have been investigated. Although our deriva¬ 
tion is the for the original quantum discord, the general Choi- 
Jamiolkowski isomorphism used here may also be applied for 
the discussion for other types of quantum discord. 

In summary, we have applied the general Choi- 
Jamiolkowski isomorphism as a convenient tool for construct¬ 
ing the transcendental equations. For the general two-qubit 
case, we have shown that the transcendental equations al¬ 
ways have a finite set of universal solutions, this result can 


be viewed as a generalization of the one get with the ARA 
algorithm. However, for some cases, the transcendental equa¬ 
tions can have solutions beside the universal ones. We also 
consider a subclass of X state, for which the transcendental 
equation may offer analytical solutions. 
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